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We generalize the classical Bass model of innovation diffusion to include a new class of agents— 
Luddites—that oppose the spread of innovation. Our model also incorporates ignorants, suscepti- 
bles, and adopters. When an ignorant and a susceptible meet, the former is converted to a susceptible 
at a given rate, while a susceptible spontaneously adopts the innovation at a constant rate. In re¬ 
sponse to the rate of adoption, an ignorant may become a Luddite and permanently reject the 
innovation. Instead of reaching complete adoption, the final state generally consists of a population 
of Luddites, ignorants, and adopters. The evolution of this system is investigated analytically and 
by stochastic simulations. We determine the stationary distribution of adopters, the time needed 
to reach the final state, and the influence of the network topology on the innovation spread. Our 
model exhibits an important dichotomy: when the rate of adoption is low, an innovation spreads 
slowly but widely; in contrast, when the adoption rate is high, the innovation spreads rapidly but 
the extent of the adoption is severely limited by Luddites. 

PACS numbers: 05.40.-a, 02.50.-r, 89.75.-k, 89.65.-s 


I. INTRODUCTION 

Models of innovation diffusion seek to understand how 
new ideas, products, or practices spread within a soci¬ 
ety through various channels [1]. Innovation may refer 
to new technologies or deviations from existing social 
norms. Rather than a single theory, innovation diffusion 
represents a theoretical framework that encompasses a 
range of social models in which the term “diffusion” can 
mean contagion, imitation, and social learning [HU. 

Many of the traditional approaches [5] to innovation 
diffusion modeling are based on a mean-field approxima¬ 
tion and are referred to as aggregate models. An influ¬ 
ential example is the seminal Bass model iHia, where 
innovation spreads as the result of either an adopter con¬ 
verting a susceptible (contagion), or through external in¬ 
fluences on susceptibles (advertising and mass media). 
The basic outcome of the Bass model is that the time de¬ 
pendence of the fraction of adopters exhibits a sigmoidal 
shape [IISHIIIS]- Thus significant adoption arises only 
after some latency period, after which complete adoption 
is quickly achieved. 

While the Bass and related models have been success¬ 
ful in fitting historic data m, there are several limita¬ 
tions of these approaches: 

• The predictive power of the Bass model is uncer¬ 
tain usms]. 

• Aggregate models are based on infinitely large, ho¬ 
mogeneous populations nans] and cannot account 
for sample-specific differences and related fluctua¬ 
tion phenomena. 

• Bass-like models do not account for behavioral pat¬ 
terns that result from social reinforcement and 
“bandwagon” pressure HTHiaEI]. 


• Aggregate models assume a “pro-innovation bias” 
and thus cannot reproduce phenomena such as in¬ 
complete adoption milZlEolET]. 

We are particularly interested in situations where in¬ 
novation can be accompanied by controversy, suspicion, 
or rejection within some social circles, potentially leading 
to incomplete adoption. As an example, mobile phones 
are owned by 90% of Americans [22] as of 2014, but their 
use is accompanied by continued health and safety con¬ 
cerns [23] • Similarly, the coverage of the measles, mumps 
and rubella vaccine in the United Kingdom reached 
92.7% in 2013-14, below the target level of 95% cover¬ 
age for herd immunity. This incomplete adoption level 
may result from doubts about vaccine effectiveness and 
safety concerns promulgated by anti-vaccination move¬ 
ments [24l [25] . Such doubts seem to persist even in the 
face of their apparently negative consequences, such as 
the measles epidemic that seemed to have its inception 
in Disneyland at the start of 2015. 

Motivated by these facts, we introduce a model 
for the diffusion of an innovation, using a statistical 
physics approach m, in which we account for the com¬ 
peting role of “Luddites” in hindering the spread of 
the innovation. Agents may either be Luddites (op¬ 
posed to innovation), “Ignorants” (no knowledge of 
the innovation), “Susceptibles” (receptive to innova¬ 
tion), or “Adopters” of the innovation. We dub this 
the LISA (Luddites/Ignorants/Susceptibles/Adopters) 
model. The main new feature of the LISA model is 
the existence of agents that reject the innovation in re¬ 
sponse to the spread of adoption. Previous work EH has 
considered the introduction of ‘resistance leaders’ who 
spread a negative response to the innovation, akin to the 
spread of a competing innovation. The LISA model dif¬ 
fers from this approach by considering agents who re¬ 
spond in particular to the rate of uptake of the inno¬ 
vation. We use the term “Luddites” in reference to the 
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19th-century Luddism movement in which English textile 
artisans protested against newly developed labor-saving 
machinery [28]. We are interested in determining how 
Luddism limits the final level of adoption and how the 
presence of Luddites leads to a trade-off between adop¬ 
tion levels and adoption times scales. 

The LISA model is defined in the next section, while 
the behavior of the model in the mean-field limit and on 
complete graphs is investigated in Section [In] Section [TV] 
focuses on the model dynamics on random graphs and 
on a one-dimensional regular lattice. For all these sub¬ 
strates, we investigate how Luddism affects the final level 
of adoption and the time scale of adoption. We also eluci¬ 
date a dichotomy between the cases of slow but relatively 
universal adoption for low values of an intrinsic innova¬ 
tion rate, and the rapid but limited spread of innovation 
that occurs in the opposite limit. Our conclusions are 
presented in Section [V| 

II. THE LISA MODEL 

As a helpful preliminary, let us review the simpler two- 
state Bass model of innovation diffusion. Here a popu¬ 
lation consists of two types of agents: susceptibles S or 
adopters A. In the Bass model, susceptibles can become 
adopters via either of two processes: 

(a) Contagion-driven conversion: a susceptible con¬ 
verts to an adopter by interacting with another 
adopter, as represented by the process S A ^ 

A “h A. 

(b) Spontaneous adoption: a susceptible converts to an 
adopter, S ^ A. 

The characteristic feature of the Bass model is that the 
adopter density exhibits a sigmoidal time dependence, 
in which the time derivative of this density has a sharp 
peak (corresponding to an inflection point in the time de¬ 
pendence of the density itself), before complete adoption 
eventually occurs [DIBITS]. 

. 

FIG. I: {Color online) Schematic depiction of the LISA 
model. An ignorant X can become a Luddite C with rate 
rA (in a mean-field setting); an ignorant can also become 
a susceptible S by contagion with rate proportional to the 
susceptible density. A susceptible spontaneously becomes an 
adopter at rate 7 . 

Our LISA model is a four-state system that consists 
of a population of N individuals that can each be in the 
states of Luddite (£), ignorant (X), susceptible (5), or 
adopter (A). Ignorant agents may either be persuaded to 
become susceptible, and thence reach the adopter state. 


or they may become a Luddite and permanently oppose 
the spread of the innovation. Specifically, the elemental 
steps of our LISA model are the following (see Fig. [^: 

(a) Contagion-driven conversion: An ignorant agent 
becomes susceptible by interacting with another 
susceptible agent. That is, X + 5 ^ S -h S with 
rate I. 

(b) Spontaneous adoption: A susceptible agent spon¬ 
taneously becomes an adopter, S ^ A with rate 7 

m- 

(c) Luddism: Ignorants may permanently reject the 
innovation and become Luddites, X ^ with a 
rate proportional to the change in the density of 
adopters in its neighborhood. 

The Luddism mechanism outlined above incorporates 
two aspects of negative behavior towards innovation. 
The first represents a fear of innovation or its conse¬ 
quences, as in the case of the historical Luddism move¬ 
ment, where the introduction of labor-saving machinery 
caused fear over job security [28|. The second is that of 
non-conformity; agents may oppose the innovation sim¬ 
ply due to its rapid increase in popularity [ 21 ] . We model 
this feature by defining the rate at which the Luddite 
density increases to be proportional to the adoption rate, 
with constant of proportionality denoted by r, the Lud¬ 
dism parameter. 

The multi-stage progression X S ^ A may also 
be viewed as a type of social reinforcement mechanism 
in which adoption follows from a succession of prompts 
from neighbors [181 IE]- The equivalent 3-state model 
with only Luddites, ignorants and adopters creates a po¬ 
larized community creates a polarized community where 
the ratio of adopters to Luddites is dependent only on the 
Luddism parameter, r. Other relevant models [27] have 
found that high levels of advertising can prompt a neg¬ 
ative response to innovation which cannot be replicated 
with only three states. The combination of a multi-stage 
progression to adoption, together with the Luddite mech¬ 
anism, arguably represents the simplest generalization of 
the Bass model that gives rise to non-trivial long-time 
state with incomplete adoption of an innovation. 


III. MEAN-FIELD DESCRIPTIONS 

We first consider the LISA model in the mean-field 
limit, where agents are perfectly mixed. The densi¬ 
ties of each type of agent are given by (X, /, S', A) = 
{NjC, Ax, A 5 , A^)/A, where Nx is the number of agents 
of type X G {£,X, S,^}, and A is the total number of 
agents. We consider the limit A ^ oc, so that all densi¬ 
ties are continuous variables and all fluctuations are neg¬ 
ligible. In this setting, the evolution of the agent densities 
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is described by the rate equations: 

L = rAl={a- 1)5'/, 

i=-{^ + ir)SI^-aSI, 

S = S{I-j), 

A = iS, 

where the dot denotes the time derivative and we define 
a = 1 + yr. Since the total density is conserved, i.e., 
// + / + 5 + 74 = 1 , the sum of these rate equations 
equals zero. A natural initial condition is a population 
that consists of a small density of susceptible agents that 
initiate the dynamics, while all other agents are ignorant; 
that is, /(O) = 1 — 5(0) = /q and 1/(0) = A(0) = 0. 

To solve these rate equations, it is useful to introduce 
the modified time variable dr = S{t)dt, which linearize 
the rate equations to 


L' = {a- 1)/, 
/' = -al, 

S'= 1-1^ 

A'=7, 


( 2 ) 


with solution 

i = —/o(l-e-“"), 
a 

I = 

° (3) 

5 = -(l-e-“") + l-Jo-7r, 

a 

A = 7T. 

There are two basic regimes of behavior that are con¬ 
trolled by the adoption rate 7 , as illustrated in Fig. 

(a) Extensive adoption. When 7 < /q, the density 
of susceptibles 5 varies non-monotonically in time 
and reaches a maximum value 5inc at an “incep¬ 
tion” time tine, after which 5 decays to 0. This 
non-monotonicity leads to a sigmoidal curve for 
the adopter density, with A increasing rapidly for 
t ^ ^inc and increasing very slowly for t > tine- The 
rescaled inception time rinc is determined by the 
criterion S' = 0, or equivalently, /(rinc) = 7 - This 
gives 

Tine = - ln(/o/7) • ( 4 ) 

a 

(b) Sparse adoption. When 7 > /q, the susceptibles 
quickly become adopters, leaving behind a substan¬ 
tial static population of ignorants and a small frac¬ 
tion of adopters, as well as Luddites. 



(a) 



(b) 


FIG. 2 : (Color online) Evolution of a realization of the LISA 
model on a complete graph of 10® nodes with Iq = 0.8 and 
Luddism parameter r = 0.9. (a) 7 = 0.3 (extensive adoption) 
(b) 7=1 (sparse adoption). Evenly distributed samples of 
the stochastic simulation (□) are indistinguishable from the 
solution of Eq. Q (solid line). The completion times for (a) 
and (b) are 60 and 17 respectively. 


We can express the densities in terms of the physical 
time t by inverting dr = S(t)dt to give t = dr'/S( t'). 
Substituting 5(r) from the third of Eqs. Q and taking 
the limits of low adoption, 7 <C 1 and a ^ 1, we have [34] 




In particular, the physical inception time tine is. 


ln(/o/7) 


Jo i 


loe-^' 


In 


lo 


(1 - Jo)7, 


( 5 ) 


(6) 


Numerical simulations of the LISA model on a large com¬ 
plete graph and numerical integration of the rate equa¬ 
tions 0 , illustrated in Fig. give results that are vir¬ 
tually indistinguishable. 


and therefore grows as ln(l/ 7 ). 

The stationary state is reached when all susceptibles 
disappear, so that no further reactions can occur. This 
gives the criterion 5(roo) = 0 which defines the value of 
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Too • By solving the third line of Eq. ^ , we obtain 

Too = -- — + -Wo ^ (7) 

'y a a \ j J 

where Wo{z) is the principal branch of the Lambert func¬ 
tion W{z)^ which is defined as the solution of 2 ; = We^. 
Here Too is a decreasing function of the adoption rate 7 , 
with Too ^ 1/7 in the high and low adoption rate regimes. 

We now determine the final densities by substituting 
Too into Eqs. Eor small adoption rate (7 ^ 1), this 
gives 

^00 = 1 - 0 ( 7 ), 

Iqo O5 

Loo^{a-l)lo = 0{y). 

Similarly, the densities at the inception time are ob¬ 
tained by substituting rinc into Eqs. (§. This yields 
A{Tinc)^S{Tinc) = Since {a-l)Io ^ 

0 ( 7 ), when 7^1 and r is finite, here the stationary 
density of adopters approximately equals the sum of the 
adopter and susceptible densities at the inception time, 
Aoo ~ ^(^inc) + 5'(rinc)- Heucc, in the low adoption rate 
regime (when r is finite), we can infer the final level of 
adoption from the adopter and susceptible densities at 
the inception time, i.e., well before the stationary state. 

The dependence of the final densities for different pa¬ 
rameter ranges is shown in Eig. Again simulation re¬ 
sults for the complete graph are indistinguishable from 
numerical integration of the rate equations. Interestingly, 
Loo varies non-monotonically on 7 when the initial state 
consists mostly of ignorants and the fixed rate of Lud¬ 
dism r is not too high, as in Eig. (top). This non¬ 
monotonic dependence on 7 can be understood by not¬ 
ing that dLooldy ^ r(l — > 0 for 7 ^ 1 and 

dLoo/dy ^ < 0 for 7 1 . We therefore ex¬ 

pect that Loo is peaked for an intermediate value of 7 on a 
range between the slow and quick adoption regimes. It is 
also worth noting that in the absence of Luddites, com¬ 
plete adoption is almost, but not completely achieved, 
since the final densities of adopters and ignorants are 
Aoo ^ 1 — loo and loo ~ see Eig.3 (bottom). 

To assess the role of finite-fluctuations on the dy¬ 
namics, we simulate the LISA model on complete graphs 
of N nodes using the Gillespie algorithm [29]. At 
long times we And that the densities of each species, 
Nx/N^ fluctuates around the corresponding mean-field 
density, with a root-mean-square fluctuation of ampli¬ 
tude ^ as expected from general properties of 

this class of reaction processes m- We also And that 
the probability distribution of Nx/N is a Gaussian of 
width of order that is centered on the mean-field 

density. We also estimate the completion time Tc for the 
system to reach its Anal state by the physical criterion 
that S{t = Tc) = 1 /A. That is, completion is defined 
by the presence of a single susceptible remaining in the 
population m- By linearizing the rate equations 0 



r 


FIG. 3: (Color online) Dependences of the final-state densi¬ 
ties Loo, loo and Aoo for a complete graph of 10 ^ nodes and 
/o = 0.9. In the top panel r = 0.9 while 7 varies, whereas 
in the bottom panel 7 = 0.3 while r varies. Simulations (□) 
in complete agreement with it with substitution 0 (solid 
line). 

around Soo = 0 , the density of susceptibles asymptoti¬ 
cally vanishes as S{t) ^ Hence, we estimate 

the mean completion time to be Tc ~ (In A )/(7 — Zoo)- 
This prediction is confirmed by our simulations. 

IV. LISA MODEL ON RANDOM GRAPHS AND 
LATTICES 

We now consider the behavior of the LISA model on 
Erdds-Renyi random graphs and one-dimensional lat¬ 
tices. We are particularly interested in uncovering dy¬ 
namics that are characterized by genuine non mean-field 
effects. 

A graph with A nodes can be represented by its A x A 
adjacency matrix A = [A^j], where Aij = 1 if nodes i and 
j are connected and 0 otherwise. We implement the LISA 
model on such a graph using the Gillespie algorithm [29] . 
The propensity for a susceptible to become an adopter 
is 7 , independent of the local environment. The propen¬ 
sity for an ignorant node i to become susceptible if it has 
Si susceptible neighbors is Si/N. The propensity of an 
ignorant node i to become a Luddite is rysijki, where 
ki = Aij is the degree (number of neighbors) of node 
i, and Si/ki is the fraction of nodes in the neighborhood 
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of i that are in the susceptible state. Thus the propen¬ 
sity of i to become a Luddite is proportional to the sum 
of its susceptible neighbors’ propensities to adopt. This 
rate encodes node i’s local knowledge of the rate of adop¬ 
tion. These reaction rates approach those of the complete 
graph, described in Section m as the average degree of 
the graph increases. 


A. Erdos-Renyi random graphs 

We first study the LISA model on the class of Erdos- 
Renyi (ER) random graphs in which an edge between 
any two nodes occurs with a fixed probability p. This 
construction leads to a binomial degree distribution for 
the ER graph in which each node has, on average, k = 
p{N — l) neighbors [30]. Under the assumption of no cor¬ 
relations between the degrees of neighboring nodes, the 
adjacency matrix may be written as Aij ^ kikj/{Nk) ~ 
k/N. The LISA dynamics on ER graphs can now be ap¬ 
proximately described by a natural generalization of the 
mean-field theory in which there are suitably defined re¬ 
action rates. In particular, if Si is the probability that a 
node i is susceptible and Ij is the probability that a node 
j is ignorant, then the density of susceptibles S evolves 
as 


Si = Si[j2iAij/N)I^-7 
j 


S[{k/N)I-j], 


since each susceptible interacts with k of its N neighbors 
on average. Thus on the ER graph there is a rescaling of 
the rate of the two-body contagion process X-\-S —5+5, 
whereas the rates of the remaining one-body processes re¬ 
main unaltered. Hence we obtain the effective rate equa¬ 
tions 


L = jrSI = (/3 - SI, 

A = jS, 

where, for later convenience, we define /3 = yr + (k/N). 

As in the case of the mean-field dynamics, the above 
equations predict two regimes of behavior (see Eig. [^: 

(a) Slow but extensive adoption (7 < klo/N). Here 
the density of 5’s peaks at a inception time tine ^ 
ln(l/ 7 ) before vanishing. 

(b) Rapid but sparse adoption (7 > klo/N). The den¬ 
sity of 5’s vanishes quickly so that the density of 
adopters and Luddites quickly reach their steady- 
state values. 




FIG. 4: {Color online) The evolution, averaged over 100 re¬ 
alizations, of the LISA model on an ER graph with N — 10^ 
nodes, k = 10 , and Iq = 0 . 8 . (a) 7 = 0 . 002 , such that 

7 < {k/N)Io and (b) 7 = 0.1 such that 7 > {k/N)Io. Shown 
are the evenly distributed samples of the stochastic simula¬ 
tion (□) and the solution of Eq. § (solid line). The Luddism 
parameter r = 0.9. 


The simulation results presented in Eig. [^indicate that 
the mean-field approximation (i correctly captures the 
main qualitative features of the dynamics on large ER 
graphs. When 7 < klo/N^ the densities of A and jC are 
characterized by a sigmoidal time dependence, whereas 
the density of 5 has a peak at the inception time tine, with 
time evolution that is slower than on complete graphs, 
since each agent has now a finite neighborhood. 


The stationary state can be determined by again not¬ 
ing that (§ becomes linear in terms of the variable 


HI 


we 


T = Jq S{t')dt'. Thus proceeding as in Section 
find the steady state by the condition S{too) = 0 . i’his 
yields 


Joo = 

l3-klN 


Lee — 




{lo - /oo) 


-^00 — 'y'^oo. 


( 9 ) 
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where now 


'^oo 


k 

N-ff3 


+ (l-/o) 


r 

'p 


^ —Wo ( 

/3 \ N'y 


( 10 ) 



FIG. 5: {Color online) Dependence of the final densities 
Loo , loo and Aoo on the average degree for ER graphs with 
iV = 10^ nodes. The simulation (□) represents an average 
over 40 model realizations for 30 randomly generated net¬ 
works. Parameters are 7 = 0.005, r = 0.9, and Iq = 0.9. The 
mean-field predictions ( 1 ^ (solid line) match the simulation 
for k >20 (see main te:?^. 


It is instructive to compare the predictions with 
the results of stochastic simulations, and also compare 
with the equivalent results for the complete graph. Fig¬ 
ure 1 ^ shows simulation results for the stationary den¬ 
sities as a function of the mean degree. These results 
confirm that the mean-field predictions correctly capture 
the functional dependence of the steady state on the vari¬ 
ous parameters. However, the mean-field predictions © 
are quantitatively accurate only when k/N is large. If 
k/N 1, the neighborhood of each agent represents a 
small fraction of the entire network, and resulting large 
demographic fluctuations invalidate the assumptions un¬ 
derlying the derivation of The dependence on 7 and 
r are qualitatively similar to those observed on complete 
graphs. 

The influence of demographic fluctuations can be 
heuristically assessed by viewing ER graphs of mean de¬ 
gree /c as a meta-population that consists of N/k patches 
each comprising a well-mixed population of size k. Ac¬ 
cording to this picture, when N ^ k ^ the number 
of agents in each of the N/k components fluctuates in 
a range of k^^‘^ about its average value. Since these are 
independent fluctuations, the noise in the whole popula¬ 
tion should have an amplitude 

which leads to fluctuations in the densities of order 
^- 1 / 2 ^ This prediction is confirmed by our simulations— 
we find that Na{oc)/N has a Gaussian probability dis¬ 
tribution around with a width that decays as 
The same behavior is observed for Loo and I^o but not 


for Soo as Soo = 0 is a requirement for the completion of 
the dynamics. 

The mean-field steady state predictions © are sum¬ 
marized in Fig. where we plot the mean-field steady 
state predictions corresponding to each pair of steady 
state densities being equal. For example, the solid curve 
corresponds to parameter values for which Aoo = ^00 • 
This carves up the ( 7 , r) parameter space into regions 
corresponding to different orderings of the steady state 
densities, labelled with Roman numerals in Fig. We 
can use these orderings to interpret, from a marketing 
perspective, whether or not these would be considered 
successful campaigns. In this context, the most desirable 
outcome would be region (I), where adopters form the 
largest steady-state group and Luddites form the small¬ 
est. Whilst adopters also form the largest group in region 
(II), Luddites form the second largest group and so this 
region could be considered a controversial success — de¬ 
spite the majority adopting, a significant number of peo¬ 
ple have responded negatively. Conversely, regions (III) 
and (IV) could both be considered controversial failures 
because Luddites form the largest groups. Regions (V) 
and (VI) represent ineffective campaigns because igno¬ 
rant s form the largest steady state groups. 



FIG. 6 : {Color online) The mean-field steady state predic¬ 
tions © over the parameter space ( 7 , r) for k/N = 0.025 and 
Iq = 0.9. The contours Loo = -^00, Loo = Aoo, and loo = Aoo 
split the domain into six regions which characterize the inno¬ 
vation (see main text). 


In summary, we have shown that the LISA dynamics 
on ER graphs can be accurately approximated by using 
mean-field assumptions, provided that the average degree 
is sufficiently high (see Fig. [^. 


B. One-dimensional lattices 

Recent controlled experiments have shown that innova¬ 
tion may spread more efficiently on clustered graphs and 
lattices than on random networks m- To understand the 
effect of regular topology on the spread of an innovation, 
and where the mean-field approximation breaks down. 
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we investigate the LISA dynamics on one-dimensional 
lattices. 

The two regimes of behavior predicted by the mean- 
field dn^r ip t ion (§ on ER random graphs (see Sec¬ 


tion 


IV A) also occur on one-dimensional lattices, despite 


the difference in topology. Specifically with k = 2 we ob¬ 
serve slow adoption for 7 < {2/N)Io and fast adoption 
for 7 > {2/N)Iq. From simulations, illustrated in Fig.7, 
we observe the following three regimes: 


(A) When 7 ^ 2/o/V, there is slow adoption as well 

as a time-scale separation. First, almost all Vs 
are converted to 5’s [32j in a time of the order of 
N^. When the lattice consists almost entirely of 
tS’s, these become adopters after a mean time of the 
order of 7 “^. Asa consequence, when 7 ^ N~^ the 
size of the adopter domains grows abruptly after a 
time of order ^ + 7 “^, when all ignorants have 

disappeared and the entire lattice is covered with 
adopters. 

(B) When 7 ^ 2/o/V, the domains of adopters grow 
initially nearly linearly in time, whereas the average 
size of I clusters remains approximately constant 
and of a comparable size to A domains. 

(C) When 7 ^ (2/V)/o, adoption occurs quickly and 
the final state is reached in a time of order 0(1/7). 
The final adopter density is limited by the forma¬ 
tion of Luddites at the ends of ignorant domains 
which prevent further conversion within each do¬ 
main. 


While the mean-field approximation predicts the 
correct regimes of behavior, the agreement is only qual¬ 
itative. In Fig. we compare typical simulations of the 
LISA model on a one-dimensional lattice with the mean- 
field predictions of ^ for the case of k = 2. The simu¬ 
lations and mean-field predictions (i systematically de¬ 
viate; the latter always overestimates the final density 
of adopters and underestimates the final density of ig¬ 
norants. This can be attributed to the topological con¬ 
straints on one-dimensional lattices. Initially the lattice 
comprises of contiguous domains of ignorants that are 
separated by domains of one or more neighboring sus- 
ceptibles. Since ignorants can only become susceptible 
if a neighbor is susceptible, domains of ignorants shrink 
at their interfaces with susceptibles. Crucially, the evo¬ 
lution of an ignorant-susceptible interface ceases if either 
the susceptible at the interface adopts or the ignorant 
at the interface becomes a Luddite. Thus in one dimen¬ 
sion both Luddites and adopters act as barriers to the 
spread of adoption, an effect that is not captured by the 
mean-field description. 

Since domains of ignorants decrease in size and evolve 
independently, we can determine analytically the ex¬ 
pected final length of ignorant domains {x) and hence 
the final fractions of each type of agent. The details 
of these calculations are given in Appendix Briefly, 




Node 


FIG. 7: {Color online) Final simulated average proportions 
of adopters (red/gray □), ignorants (green/dark gray □) and 
Luddites (blue/black □) for varying values of 7 , averaged over 
100 simulations. Theoretical predictions using ignorant do¬ 
main length (see Appendix [A| for details) are overlaid (solid 
line). Parameters are N = 1000, r = 0.5. Initially igno¬ 
rants and susceptibles are randomly distributed, with densi¬ 
ties Iq = 0.8 and So = 0.2. The three regimes discussed in the 
text are separated by dashed lines corresponding to regions 
where (2/A)/o ^ 7 and {2/N)Io ^ 'y. Typical realizations of 
the model for A = 100 in each of the three regimes are given 
(bottom). On the vertical axis the iteration corresponds to a 
single step of the Gillespie algorithm, with one reaction taking 
place per iteration. 


we first determine the probability Pn{'m) that a domain 
of ignorants of initial length n shrinks by m. We then 
use Pn{Tn) to calculate the expected final length of ig¬ 
norant domains {x) and the final fraction of ignorants. 
Since Luddites only form at the boundaries of ignorant 
domains, we are able also to determine the expected fi¬ 
nal fraction of Luddites and hence, using the conserva¬ 
tion relation L I S A = 1 , the final fraction of 
adopters. The resulting final fractions of each type of 
agent are plotted in Fig. and agree extremely well with 
the corresponding numerical simulations. In principle, 
this method allows us to derive explicit formulas for the 
final fractions of each agent; however, in practice these 


















FIG. 8 : {Color online) Time dependence of the densities in 
each state for a one-dimensional lattice of size iV = 10^ av¬ 
eraged over 100 realizations. The corresponding mean-held 
predictions from Eq. with k — 2 (solid line) deviate dra¬ 
matically from the simulation samples (□). The parameters 
are 7 = 0.005, r = 0.9, and Iq = 0.8. 


formulas prove cumbersome. 


V. DISCUSSION CONCLUSION 

luuovatious are ofteu accompauied by societal debates 
aud coutroversies that may lead to divisious betweeu 
adopters of au iuuovatiou aud those who permaueutly 
reject that iuuovatiou. Cousequeutly, iuuovatious are 
rarely adopted by the whole populatiou, as various exam¬ 
ples, raugiug from techuology to mediciue, demoustrate. 
Classical models of iuuovatiou diffusiou, such as that pro¬ 
posed by Bass, assume a “pro-iuuovatiou bias” aud pre¬ 
dict the complete adoptiou of iuuovatious. 

Motivated by these cousideratious we have iutroduced 
a multi-stage geueralizatiou of the Bass model, the LISA 
model, that does uot uuavoidably lead to complete adop¬ 
tiou. The maiu uew feature of our model is the iutroduc- 
tiou of Luddites that permaueutly oppose the spread of 
iuuovatiou iu their ueighborhood. lu the LISA model, ig- 
uoraut iudividuals cau successively become susceptibles 
aud theu adopters, or turu to Luddism iu respouse to a 
high rate of adoptiou aud permaueutly reject the iuuo¬ 
vatiou. 

We carried out a detailed aualysis of the properties 
of the LISA model ou complete graphs aud Erdos-Reuyi 
raudom graphs, as well as ou oue-dimeusioual lattices. Iu 
particular, we focused ou the steady states aud comple- 
tiou time (time to reach statiouarity). We showed that 
siguificaut iusights cau be gaiued from a simple meau- 
field aualysis that aptly captures the qualitative aspects 
of the two basic regimes of the LISA dyuamics. Wheu 
the rate of adoptiou is low, the populatiou slowly cou- 
verges to a dual state that cousists of a high couceutra- 
tiou of adopters. Iu the couverse case, the statiouary 
state is reached much more quickly, but the dual fractiou 


of adopters is much lower aud is severely limited by the 
siguidcaut deusities of Luddites aud iguorauts. 

Siuce most models of iuuovatiou didusiou are formu¬ 
lated at meau-deld level, au importaut aspect of this work 
has also beeu to reveal the limitatious of the meau-deld 
approximatiou. Iu particular, for Erdos-Reuyi raudom 
graphs with low meau degree aud oue-dimeusioual lat¬ 
tices, the meau-deld approximatiou proves iuaccurate. 
This is due to the formatiou of Luddites which isolate 
domaius of iguorauts from the iuuovatiou, au edect par¬ 
ticularly appareutly iu oue dimeusiou. It would be worth¬ 
while to iuvestigate the LISA model ou modular uet- 
works, where Luddism has the poteutial to block the 
spread of iuuovatiou to eutire commuuities. Iu addi- 
tiou to the work described iu this paper, we also fouud 
that the meau-deld approximatiou proves better ou two- 
dimeusioual lattices thau ou oue-dimeusioual lattices. 

Iu summary, the LISA model is a simple, but uou- 
trivial, iuuovatiou diffusiou model that accouuts for the 
possibility that the promotiou of au iuuovatiou may be 
tempered by the alieuatiou of some iudividuals. These 
iu turu affect the spread of the iuuovatiou. luterestiugly, 
our model outliues two possible marketiug sceuarios: If 
oue is iuterested iu reachiug a high level of adoptiou theu 
this cau ouly be achieved over loug time scales, siuce the 
rate of adoptiou must be low. However, if the priority is 
to attaiu a duite level of adoptiou as quickly as possible 
regardless of the alieuatiou that this may cause, theu a 
high rate of adoptiou is preferable. 


VI. ACKNOWLEDGEMENTS 

This work is supported by au EPSRC ludustrial Case 
Studeutship Graut No. EP/L50550X/1. SR is supported 
iu part by NSF Graut No. DMR-1205797. Partial fuud- 
iug from Bloom Ageucy iu Leeds U.K. is also gratefully 
ackuowledged. 


Appendix A: Analysis of one-dimensional dynamics 

Iu this appeudix we describe the calculatiou of the d- 
ual fractious of each type of ageut ou oue-dimeusioual 
lattices. These results are compared with simulatious iu 
Fig. of Sectiou 


IV B 


1. Analysis of ignorant domains 

luitially, the uodes ou the oue-dimeusioual lattice are 
either iguoraut, with probability Jq, or susceptible, with 
probability So = 1 — Jq. Thus the iuitial coudguratiou 
cousists of couuected domaius of iguoraut uodes bordered 
by susceptibles. Moreover, siuce iguorauts cau ouly be¬ 
come susceptible if a ueighbor is susceptible, domaius of 
iguorauts ouly evolve at their iguoraut-susceptible iuter- 
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faces. We will refer to these as “active interfaces”. At an 
active interface one of three events can occur: 

• The ignorant node becomes susceptible, thus reduc¬ 
ing the domain length by one, with probability 

l/N 

“ l/A^ + r 7/2 + 7 ' 

• The ignorant node becomes a Luddite, thus reduc¬ 
ing the length of the domain by one and causing 
the interface to become inactive, with probability 

r 7/2 

1 /A^ + r 7 / 2 - 1-7 * 

• The susceptible node becomes an adopter, thereby 
terminating the interface evolution, with probabil¬ 
ity 

= 7 

l/A^ + r 7/2 + 7 ’ 

For an isolated ignorant node with two susceptible 
neighbors, these probabilities respectively become 

. _ ‘2./N 

“ 2 /A^ + r 7 + 7 ’ 

'^7 

2 /A^ + r7-h7’ 

. ^ 7 

2 /A^-l-r 7 + 7 

Let Qriijy^) be the probability that a domain of igno- 
rants of initial length n with a single ignorant-susceptible 
interface has a final length n — with 0 < m < n. We 
can determine Qn(^) as follows: If the final length of 
ignorants is n — m, with 0 < m < n, then either m igno¬ 
rant nodes must become susceptible before a susceptible 
node at the interface adopts, or m — 1 ignorant nodes 
must become susceptible before an ignorant node at the 
interface becomes a Luddite. These events occur with 
probabilities paP^ and PlP^~^ respectively. Using sim¬ 
ilar reasoning for the cases m = 0 and m = n, we thus 
find 

{ Pa if m = 0 

pap's +PlP^~^ if 0 < m < n . (Al) 
Ps+PlpT^ if m = n 

By summing over m, it can be shown that Qn(^) is nor¬ 
malized. 

We now consider the case where a connected region 
of n ignorant nodes initially has two ignorant-susceptible 
interfaces. The probability Pn{Tn) that a region of ig¬ 
norants of initial length n with two active interfaces has 
final length n — m is given by the recursion relation 

Pn{m) = Qn{m)pA + Qn-i{m - 1 )pl 

+ Pn-i{m-l)ps, (A 2 ) 


where the terms Qra{m) are given by Equation ( |A 2 [ ) 

captures the three possible events that can occur at the 
interface. If a susceptible node at the interface adopts, 
which occurs with probability pa^ then the region of ig¬ 
norants only has one remaining active interface left and 
there will be n — m remaining ignorants with probability 
(5n(^): s-s given in ( |A 1 [ ). If an ignorant node at the in¬ 
terface becomes a Luddite, which occurs with probability 
Pl, then again the region of ignorants will only have one 
active interface. Since there will be one ignorant less the 
probability there will be n — m remaining ignorants is 
(5n-i(^ ~ !)• Finally, if an ignorant node at the bound¬ 
ary becomes susceptible, which occurs with probability 
ps^ then the probability that there are n — m ignorants 
remaining is the same as if we had started with n — 1 

ignorant nodes, i.e. — 1 ). _ 

To solve the recursion relation ( |A 2 [ ) we need Pn{0) and 
Pi (1). The probability that a region of ignorants of initial 
length n remains of length n is given by 


Pn{ 0 ) = 


Papa if n = 1 
Pa 


^ if n > 1 


Also, the probability that a single ignorant node that 
initially has two susceptible neighbors becomes a suscep¬ 
tible or Luddite is given by 

Pi(l) = pa{pl + Ps) +PlP PS- 


Thus the solution to the recursion relation (A 2 ) for 0 < 
m < n — 1 is given by 

P„(m) = (m + l)p\p's + ‘2mpAPLP^~^ 

+ (m - 1)pIp^~^. 

For m = n — 1 we have 

P„(n -l)=PA\pA + {n- 1 )pa] Ps~^ 

+ 2 (n - 1)paPlpT'^ + (n - 2)plp^~^, 
and for m = n we have 


Pn{n) = \Pa{pl+Ps)+Pl+Ps]Ps ^ 

+ (n - 1) (pAPs + ‘^PlpT^ + pIpT"^) • 

Again it is possible to check, by summing ( |A 2 | ) over m 
and solving the resulting recursion relation, that Pn{Tn) 
is normalized. 

We can use Pn{'m) to calculate the expected final 
length of ignorant domains {x). First note that since Jq is 
the initial probability of being ignorant, the probability 
that a domain of ignorants initially has length n > 0 is 
given by po{n) = Iq~^So for large N. Thus we find that 

N N n 

(^) = - '^po{n)'^lPn{l)- 

n=0 n=0 1=0 

In principle, we may use the above to obtain an explicit 
expressio n fo r {x). In practice, however, we use the solu¬ 
tions to ( |A 2 [ ) to calculate {x) numerically. 
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2. Calculation of population densities 

Initially, the mean number of ignorants is given by IqN 
and so dividing by the mean length of ignorant domains, 
1/(1 — /o), yields the expected number of ignorant do¬ 
mains, (1 — /o)/oA^. Thus the final density of ignorants 
is 

^CX) = (1 — ^o)^o(^)- 

The probability that an ignorant domain survives is 

CX3 

q = l- '^po{n)Pn{n). 

n=0 


in Luddites, it is possible to determine that the expected 
number of Luddites that arise when a domain of initial 
size n > 1 vanishes is given by 

In = {paPl + Pl)P s~^ + (n - 1) {2pLPs~^ +p\pT'^) 

and /i = paPl Thus the expected number of Lud¬ 

dites that arise from domains of ignorants that vanish 
is 

oo 

?7o = (1 - Io)IoN'^pQ{n)ln. (A4) 

n=0 


Surviving ignorant domains have two interfaces, which 
are either ignorant-adopter or ignorant-Luddite, with 
probabilities pa/{pl ) and Pl/{pl-\-PA ), respectively. 
Thus the expected number of Luddites at the interfaces 
of non-vanishing ignorant domains is given by 

r]+ = - Io)IoN. (A3) 

PL ^ PA 

It is also possible for Luddites to arise when a domain 
vanishes. By identifying the terms in Pn{ri) that result 


Summing Eqs. (A3) and (A4) and dividing by N we ar¬ 
rive at the final density of Luddites 


Lo,=/o(l-/o) q+y^Po{n)ln 

\PL+PA 

Since the dynamics cease when S' = 0, the number of 
adopters can be found using the conservation law Aoo = 

1 Tqo IqO' 



[1] E. M. Rogers, Diffusion of Innovations (Free Press, New 
York, 2003). 

[2] D. L. Kincaid, J. of Health Comm. 9, 37 (2004). 

[3] B. Ryan and N. Gross, Rural Sociology 8, 15 (1943). 

[4] J. Coleman, E. Katz, and H. Menzel, Sociometry 20, 253 
(1957). 

[5] L. A. Fourt and J. W. Woodlock, Journal of Marketing 
25 31 (1960); E. Mansfield, Econometrica 29, 741 (1961). 

[6] F. M. Bass, Management Science 15, 215 (1969). 

[7] F. M. Bass, J. Business 53, S51 (1980). 

[8] V. Mahajan, E. Muller, and F. M. Bass, Journal of Mar¬ 
keting 54, 1 (1990). 

[9] F. M. Bass, Management Science 50, 1833 (2004). 

[10] W. J. Hopp, Management Science 50, 1763 (2004). 

[11] M. G. Dekimpe, P. M. Parker and M. Sarvary, Technol. 
Forecast. Soc. Change 57, 105 (1998); S. Sundqvista, 
L. Franka and K. Puumalainen, J. Bus. Res. 58, 107 
(2005); C. Michalakelis, D. Varoutas and T. Sphicopou- 
los, Telecommun. Policy 32, 234 (2008): J. Lim, C. Nam, 
S. Kim, H. Rhee, E. Lee and H. Lee, Telecommun. Policy 
36, 858 (2012); J. L. Toole, M. Cha, and M. C. Gonzalez, 
PLoS One 7, e29528 (2012). 

[12] E. Kiesling et al., CEJOR 20 , 183 (2012). 

[13] H. Peyton Young, American Economic Review, 99, 1899 
(2009). 

[14] F. Sultan, J. U. Farley, D. R. Lehmann DR, Journal of 
Marketing 27 70 (1990). 

[15] R. M. Heeler, T. P. Hustad, Management Science 26, 
1007 (1980); C. Van den Bulte, G. L. Lilien, Marketing 
Science 16, 338 (1997); R. Kohli, D. R. Lehmann, J. Pae, 
J. Product Innovation Management 16, 134 (1999). 


[16] M Hohnisch, S. Pittnauer, and D. Stauffer, Industrial 
and Corporate Change 17, 1001 (2008). 

[17] D. Centola, Science 329, 1194 (2010). 

[18] D. Centola, R. Wilker, and M. W. Macy, Am. J. Sociol. 
110, 1009 (2005); D. Centola, V. M Eguiluz, and M. W. 
Macy, Physica A 374, 449 (2007). 

[19] P. S. Dodds and D. J. Watts, Phys. Rev. Lett. 92, 
218701 (2004); H. P. Young, Amer. Econ. Rev. 99, 1899 
(2009); P. L. Krapivsky, S. Redner, and D. Volovik, 

J. Stat. Mech. P12003 (2011). 

[20] J. Goldenberg, B. Libai, S. Solomon, N. Jan N, and D. 
Stauffer, Physica A 284, 335 (2000); D. Strang D and 
M. W. Macy, American Journal of Sociology 107, 147 
( 2001 ). 

[21] P. S. Tolbert and L. G. Zucker, Admin. Science Quar¬ 
terly 28, 22 (1983); E. Abrahamson and L. Rosenkopf, 
Acad. Management Review 18, 487 (1993); L. Rosenkopf 
and E. Abrahamson, Comp, and Math. Organization 
Theory 5, 361 (1999). 

[22] http://www.pewinternet.org/data-trend/mobile/device- 
ownership / 

[23] http://interphone.iarc.fr/UICC_Report_FinaL03102011.pdf 

[24] R. M. Wolfe and L. K. Sharp, BMJ 325, 430 (2002). 

[25] http://www.hscic.gov.uk/catalogue/PUB 14949/nhs- 
immu-stat-eng-2013-14-rep.pdf 

[26] D. J. Watts and P. S. Dodds, J Cons Res 34, 441 (2007); 

G. Kocsis and F. Kun, J. Stat. Mech. P10014 (2008); 

M. Lin, N. Li, Physica A 389, 473 (2010); G. Pegoretti, 

F. Rentocchini, G. V. Marzetti, J. Econ. Interact. Co¬ 
ord. 7, 145 (2012); N. J. McCullen, A. M. Rucklidge, 

C. S. E. Bale, T. J. Foxon, and W. F. Gale, SIAM 







11 


J. Applied Dynamical Systems 12, 515 (2013); S. Mel¬ 
nik, J. A. Ward, J. P. Gleeson, M. A. Porter, Chaos, 

23, 013124 (2013); K. Sznajd-Weron, J. Szwabinski, R. 
Weron, T. Weron, J.Stat.Mech. (2014) P03007; P. Przy- 
byla, K. Sznajd-Weron, R. Weron, Advances in Complex 
Systems 17 (2014) 1450004 

[27] S. Moldovan and J. Goldenberg, Technological Forecast¬ 
ing & Social Change 71 (2004) 425442. 

[28] http://www.nationalarchives.gov.uk/education/politics/g3/ 

[29] D. T. Gillespie, J. Comput. Phys. 22, 403 (1976). 

[30] M. Newman, Networks: An Introduction (Oxford Uni¬ 
versity Press, Oxford, 2010); M. O. Jackson, Social and 
Economic Networks (Princeton University Press, Prince¬ 
ton, 2010). 

[31] N. G. Van Kampen, Stochastic Processes in Physics 
and Chemistry (Elsevier, 2007); C. Gardiner, Stochastic 
Methods (Springer, 2010). 


[32] Nonequilibrium Statistical Mechanics in One Dimension^ 
edited by V. Privman (Cambridge University Press, New 
York, 1997); S. Redner, A Guide to First-Passage Pro¬ 
cesses (Cambridge University Press, New York, 2001); 
P. L. Krapivsky, S. Redner and E. Ben-Naim, A kinetic 
view of statistical physics (Cambridge University Press, 
New York, 2010). 

[33] Adoption could also occur by contagion, according to 
S -\- A ^ A-\- A. Yet, this two-body process would yield 
similar features as our LISA model, but would be tech¬ 
nically more tedious to handle. 

[34] Here the term —yr' has been neglected. This approxi¬ 
mation is legitimate since t is integrated from 0 to r 
Too ~ 1/7 and therefore yr' ^ 1 in the regime being con¬ 
sidered. A similar reasoning, with yr < 7ln(/o/7) ^ 1, 
leads to § when 7 ^ 1 . 


